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Within a 2D model of Josephson junction arrays (created by 2D network of twin boundary dislocations 
with strain fields acting as an insulating barrier between hole- rich domains in underdoped crystals), a few novel 
effects expected to occur in intrinsically granular material are predicted including: (i) Josephson chemomag- 
netism (chemically induced magnetic moment in zero applied magnetic field) and its infiuence on a low-field 
magnetization (chemically induced paramagnetic Meissner effect), and (ii) magnetoconcentration effect (cre- 
ation of oxygen vacancies in applied magnetic field) and its infiuence on a high-field magnetization (chemically 
induced analog of "fishtail" anomaly). The conditions under which these effects can be experimentally mea- 
sured in non-stoichiometric high- Tc superconductors are discussed. 

PACS: 61.72.-y, 74.25.Ha, 74.50.-fr, 74.72.-h 
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1. Introduction. Recent imaging of the granu- 
lar structure in underdoped Bi2Sr2CaCu20s+s crys- 
tals [Q, revealed an apparent segregation of its elec- 
tronic structure into superconducting domains (of the 
order of a few nanometers) located in an electronically 
distinct background. In particular, it was found that at 
low levels of hole doping {S < 0.2), the holes become 
concentrated at certain hole- rich domains. (In this re- 
gard, it is interesting to mention a somewhat similar 
phenomenon of " chemical localization" that takes place 
in materials, composed of atoms of only metallic ele- 
ments, exhibiting metal-insulator transitions Q.) Tun- 
neling between such domains leads to intrinsic granu- 
lar superconductivity (GS) in high- Tc superconductors 
(HTS). Probably one of the first examples of GS was 
observed in Y Ba2Cu^0'j^s single crystals in the form 
of the so-called "fishtail" anomaly of magnetization 
The granular behavior has been related to the 2D clus- 
ters of oxygen defects forming twin boundaries (TBs) or 
dislocation walls within CuO plane that restrict super- 
current flow and allow excess flux to enter the crystal. 
Indeed, there are serious arguments to consider the TB 
in HTS as insulating regions of the Josephson SIS-type 
structure. An average distance between boundaries is 
essentially less than the grain size. In particular, the 
networks of localized grain boundary dislocations with 
the spacing ranged from lOnm to lOOnm have been ob- 
served 1^ which produce effectively continuous normal 
or insulating barriers at the grain boundaries. It was 
also verified that the processes of the oxygen ordering 
in HTS leads to the continuous change of the lattice pe- 
riod along TB with the change of the oxygen content. 



Besides, a destruction of bulk superconductivity in these 
non-stoichiometric materials with increasing the oxygen 
deficiency parameter 5 was found to follow a classical 
percolation theory ji) . 

In addition to their importance for understanding 
the underlying microscopic mechanisms governing HTS 
materials, the above experiments can provide rather ver- 
satile tools for designing chemically-controlled atomic 
scale Josephson junctions (JJs) and their arrays (JJAs) 
with pre-selected properties needed for manufacturing 
the modern quantum devices H, |^. Moreover, as 
we shall see below, GS based phenomena can shed 
some light on the origin and evolution of the so-called 
paramagnetic Meissner effect (PME) which manifests 
itself both in high- Tc and conventional superconduc- 
tors 1^, P and is usually associated with the presence 
of TT-junctions and/or unconventional (d-wave) pairing 
symmetry. Since recently, much attention has been paid 
to both experimental and theoretical study of PME re- 
lated effects using specially designed S'F5'-type junc- 
tions §, 0. 

In this Letter, within a 2D model of JJAs (cre- 
ated by a regular 2D network of TB dislocations), we 
discuss a possibility of a few novel interesting effects 
which are expected to occur in intrinsically granular 
non-stoichiometric material. In particular, we shall 
consider (i) Josephson chemomagnetism (chemically in- 
duced magnetic moment in zero applied magnetic field) 
and its influence on a low-field magnetization (chemi- 
cally induced PME), and (ii) magnetoconcentration ef- 
fect (creation of extra oxygen vacancies in applied mag- 
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netic field) and its influence on a high-field magnetiza- 
tion (chemically induced analog of " fishtail" anomaly) . 

2. The scenario. As is well-known, the presence of 
a homogeneous chemical potential /i through a single J J 
leads to the AC Josephson effect with time dependent 
phase difference dcji/dt = fi/h. In this paper, we will 
consider some effects in dislocation induced JJ caused 
by a local variation of excess hole concentration c(x) 
under the chemical pressure (described by inhomoge- 
neous chemical potential /x(x)) equivalent to presence 
of the strain field of 2D dislocation array e(x) forming 
this Josephson contact. 

To understand how GS manifests itself in non- 
stoichiometric crystals, let us invoke an analogy with 
the previously discussed dislocation models of grain- 
boundary Josephson junctions (GBJJs) (see, e.g., 
[l2| and further references therein). Recall that under 
plastic deformation, grain boundaries (GBs) (which are 
the natural sources of weak links in HTS), move rather 
rapidly via the movement of the grain boundary dislo- 
cations (GBDs) comprising these GBs. Using the above 
evidence, in the previous paper we studied nu- 
merous piezomagnetic effects in granular superconduc- 
tors under mechanical loading. At the same time, ob- 
served |l| in HTS single crystals regular 2D 
dislocation networks of oxygen depleted regions (gener- 
ated by the dissociation of < 110 > twinning disloca- 
tions) with the size do of a few Burgers vectors, forming 
a triangular lattice with a spacing d> do ranging from 
lOnm to lOOnm, can provide quite a realistic possibil- 
ity for existence of 2D Josephson network within CuO 
plane. Recall furthermore that in a d-wave orthorhom- 
bic YBCO crystal TBs are represented by tetragonal 
regions (in which all dislocations are equally spaced by 
do and have the same Burgers vector a parallel to y- 
axis within CuO plane) which produce screened strain 
fields 



e(x) = eoe^l^'l/* with | x | = 



r 



Though in Y Ba2Cu30'/^s the ordinary oxygen dif- 
fusion D = Doe^^'^/^^'^ is extremely slow even near 
Tc (due to a rather high value of the activation energy 
Ud in these materials, typically Ud — leV), in under- 
doped crystals (with oxygen-induced dislocations) there 
is a real possibility to facilitate oxygen transport via the 
so-called osmotic (pumping) mechanism |16, |l^ which 



relates a local value of the chemical potential (chemical 
pressure) /i(x) = /i(0) + V/i • x with a local concen- 
tration of point defects as follows c(x) = g-M(x)/fcBT^ 
Indeed, when in such a crystal there exists a nonequi- 
librium concentration of vacancies, dislocation is moved 
for atomic distance a by adding excess vacancies to the 
extraplane edge. The produced work is simply equal 
to the chemical potential of added vacancies. What is 



important, this mechanism allows us to explicitly in- 
corporate the oxygen deficiency parameter S into our 
model by relating it to the excess oxygen concentration 
of vacancies c„ = c(0) as follows 6 = 1 — Cy. As a re- 
sult, the chemical potential of the single vacancy reads 
/ij, = /i(0) = —ksT \og{l — (5) ~ ksTS. Remarkably, 
the same osmotic mechanism was used by Gurevich 
and Pashitskii to discuss the modification of oxy- 
gen vacancies concentration in the presence of the TB 
strain field. In particular, they argue that the change 
of e(x) under an applied or chemically induced pres- 
sure results in a significant oxygen redistribution pro- 
ducing a highly inhomogeneous filamentary structure 
of oxygen-deficient nonsuperconducting regions along 
GB (for underdoped superconductors, the vacan- 
cies tend to concentrate in the regions of compressed 
material). Hence, assuming the following connection 
between the variation of mechanical and chemical prop- 
erties of planar defects, namely /i(x) = Kilo^i^) (where 
rio is an effective atomic volume of the vacancy and K 
is the bulk elastic modulus), we can study the proper- 
ties of TB induced J Js under intrinsic chemical pressure 
V/x (created by the variation of the oxygen doping pa- 
rameter 6). More specifically, a single SIS type junc- 
tion (comprising a Josephson network) is formed around 
TB due to a local depression of the superconducting or- 
der parameter A(x) oc e(x) over distance do produc- 
ing thus a weak link with (oxygen deficiency S depen- 
dent) Josephson coupling J{6) — e(x)Jo — Jo{6)e^^^^^'^° 
where JoiS) = cqJo = {fJ'v/K^Q)Jo (here Jq oc Aq/Ru 
with Rn being a resistance of the junction). Thus, the 
considered here model indeed describes chemically in- 
duced GS in underdoped systems (with S ^ 0) because, 
in accordance with the observations, for stoichiometric 
situation (when (5 ~ 0), the Josephson coupling J{d) ~ 
and the system loses its explicitly granular signature. 

3. The model. To adequately describe chemomag- 
netic properties of an intrinsically granular supercon- 
ductor, we employ a model of 2D overdamped Joseph- 
son junction array which is based on the well known 
tunneling Hamiltonian 



N 



7^(0 =^J^j[l- COS 0,j(t)] 



(1) 



and introduces a short-range (nearest-neighbor) inter- 
action between N junctions (which are formed around 
oxygen-rich superconducting areas with phases 0i(t)), 
arranged in a two-dimensional (2D) lattice with co- 
ordinates Xi = {xi,yi). The areas are separated 
by oxygen-poor insulating boundaries (created by TB 
strain fields e(xij)) producing a short-range Josephson 
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coupling Jij = Jo{6)e~^^'^^^^. Thus, typically for granu- 
lar superconductors, the Josephson energy of the array 
varies exponentially with the distance = — Xj be- 
tween neighboring junctions (with d being an average 
junction size). 

If, in addition to the chemical pressure V/i(x) = 
XiloVe(x), the network of superconducting grains is 
under the influence of an applied frustrating magnetic 
field B, the total phase difference through the contact 
reads 



,W = < + -;-(xy An, 



B 



V/i • Xiji 



(2) 



where is the initial phase difference (see below). 



with X,- 



(x; -|- Xj)/2, and w 



2Xl(T) + 1 with Xl being the London penetration depth 
of superconducting area and / an insulator thickness 
(which, within the discussed here scenario, is simply 
equal to the TB thickness 0). 

To neglect the influence of the self-fleld effects in a 
real material, the corresponding Josephson penetration 
length Xj — yj $0 j^i^ix^jcW must be larger than the 
junction size d. Here jc is the critical current density 
of superconducting (hole-rich) area. As we shall see be- 
low, this condition is rather well satisfied for HTS single 
crystals. 

4. Chemomagnetism. Within our scenario, the 
sheet magnetization M of 2D granular superconductor 
is defined via the average Josephson energy of the array 



<n>= 



dt 



-H(x, t) 



as follows 



M(B,<5) = - 



d<n> 



(3) 



(4) 



where s = 2Trd^ is properly defined normalization area, 
r is a characteristic Josephson time, and we made a 
usual substitution J^ij ^ « / d'^xA(x,t) valid 

in the long- wavelength approximation [l8| . 

To capture the very essence of the superconducting 
analog of the chcmomagnetic effect, in what follows we 
assume for simplicity that a stoichiometric sample (with 
5 ~ 0) does not possess any spontaneous magnetization 
at zero magnetic field (that is M(0, 0) = 0) and that its 
Meissner response to a small applied field B is purely 
diamagnetic (that is M{B,0) ~ —B). According to 
Eq.(4), this condition implies = 27rm for the initial 
phase difference with m = 0,±l,±2,... 

Taking the applied magnetic field along the c-axis 
(and normal to the CuO plane), that is B = (0,0, B), 
we obtain finally 

^'^(^'^) = -"^"(^) (i + 6^)('i;V6.n 




Fig.0. The magnetization M{B,5)/Mo{0) as a function 
of applied magnetic field B/Bo, according to Eq.(5), for 
different values of oxygen deficiency parameter: 5 ~ 
(solid line), S — 0.05 (dashed line), and 5 = 0.1 (dotted 
line). Inset: 5 induced magnetization A/(0, 5)/Mo(0) in 
a zero applied magnetic field (chemomagnetism). 



for the chemically-induced sheet magnetization of the 
2D Josephson network. 

Here Mo{S) = Jo{5)/Bq with Jo{5) defined earlier 
(in what follows, Mo(0) is Mo{5 ~ 0)), b = B / Bo, and 

= B^/Bn ~ {kBTT/h)S where B^((5) = {fivT/h)Bn 
is the chemically-induced contribution (which disap- 
pears in optimally doped systems with 5 ~ 0), and 
Bq — ^o/wd is a characteristic Josephson field. 

Fig. shows changes of the initial (stoichiometric) 
diamagnetic magnetization M/Mq (solid line) with oxy- 
gen deficiency 5. As is seen, even relatively small val- 
ues of S parameter render a low field Meissner phase 
strongly paramagnetic (dotted and dashed lines). The 
inset of Fig. |l| presents a true chcmomagnetic effect 
with concentration (deficiency) induced magnetization 
M{0,6) in zero magnetic field. According to Eq.(5), 
the initially diamagnetic Meissner effect turns paramag- 
netic as soon as the chemomagnetic contribution B^{d) 
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Fig.g. Magnetic field dependence of tfie oxygen defi- 
ciency parameter 5{B) (magnetoconcentration effect). 

exceeds an applied magnetic field B. To see whether 
this can actually happen in a real material, let us esti- 
mate a magnitude of the chemomagnetic field Bf^ . Typi- 
cally I, H], for HTS single crystals Al(0) « ISOnm and 
d lOnm, leading to Bq ~ 0.5T. Using t ~ h/ fiy and 
jc = 10^"yl/m^ as a pertinent characteristic time and 
the typical value of the critical current density, respec- 
tively, we arrive at the following estimate of the chemo- 
magnetic field Bf,{S) ~ 0.5Bq for 5 = 0.05. Thus, the 
predicted chemically induced PME should be observable 
for applied magnetic fields B ~ 0.5i?o — 0.25T (which 
are actually much higher than the fields needed to ob- 
serve the previously discussed piezomagnetism and 
stress induced PME in high-Tc ceramics). Notice that 
for the above set of parameters, the Josephson length 
A J ~ 1/iTO, which means that the assumed in this paper 
small-junction approximation (with d <^ Xj) is valid 
and the so-called "self- field" effects can be safely ne- 
glected. 

5. Magnetoconcentration effect. So far, we ne- 
glected a possible field dependence of the chemical po- 
tential fj,y of oxygen vacancies. However, in high enough 
applied magnetic fields B, the field-induced change 
of the chemical potential A/x„(i3) = fiv{B) — /Li„(0) 



becomes tangible and should be taken into account. 
As is well-known 0, in a superconducting state 
A^iy{B) = ~M{B)B/n, where M{B) is the corre- 
sponding magnetization, and n is the relevant carriers 
number density. At the same time, within our sce- 
nario, the chemical potential of a single oxygen va- 
cancy /ii, depends on the concentration of oxygen va- 
cancies (through deficiency parameter S). As a result, 
two different effects are possible related respectively to 
magnetic field dependence of Hv{B) and to its depen- 
dence on magnetization ^y{M). The former is nothing 
else but a superconducting analog of the so-called mag- 
netoconcentration effect (which was predicted and ob- 
served in inhomogeneously doped semiconductors pl[|) 
with field-induced creation of oxygen vacancies Cv{B) — 
Cy{0)exp{—Afiy{B)/kBT), while the latter (as we shall 
see in the next Section) results in a "fishtail" -like be- 
havior of the magnetization. Let us start with the mag- 
netoconcentration effect. Figure ^depicts the predicted 
field-induced creation of oxygen vacancies 5{B) — 1 — 
Cy{B) using the above-obtained magnetization M {B, S) 
(see Fig. Q and Eq.(5)). We also assumed, for simphc- 
ity, a complete stoichiometry of the system in a zero 
magnetic field (with S{0) — I — Cu(0) = 0). Notice that 
5{B) exhibits a maximum at Sc — 0.23 for applied fields 
B = Bq (in agreement with the classical percolative 
behavior observed in non-stoichiometric YBa2C'u30'/-^s 
samples ^, ^ |l5|). Finally, let us show that in under- 
doped crystals the above-discussed osmotic mechanism 
of oxygen transport is indeed much more effective than 
a traditional diffusion. Using typical YBCO parame- 
ters [Q, eo = 0.01, flo = Oq with oq — 0.2nm, and 
K — 115GPa, we have fiy{0) — eoKflo ~ ImeV for a 
zero-field value of the chemical potential in HTS crys- 
tals, which leads to creation of excess vacancies with 
concentration c„ (0) = e^^"(°)/'''s^ ~ 0.75 (equivalent to 
a deficiency value of (5(0) ~ 0.25) at T = Tc, while the 
probability of oxygen diffusion in these materials (gov- 
erned by a rather high activation energy Ud — leV) 
is extremely slow under the same conditions because 
D oc 6"^'*/'^^-^'= <^C 1. On the other hand, the change 
of the chemical potential in applied magnetic field can 
reach as much as |2§ A/i„(B) ~ O.bmeV for B = 0.5T, 
which is quite comparable with the above-mentioned 
zero-field value of fiy (0) . 

6. Analog of "fishtail" anomaly. Let us turn 
now to the second effect related to the magnetization 
dependence of the chemical potential fiy{M{B)). In 
this case, in view of Eq.(2), the phase difference will 
acquire an extra M{B) dependent contribution and as 
a result the r.h.s. of Eq.(5) will become a nonlin- 
ear functional of M{B). The numerical solution of 
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B/B„ 

Fig.^. A "fishtail" -like behavior of magnetization mf = 
M{B,5{B))/Mo{0) in applied magnetic field B/Bo in 
the presence of magnetoconcentration effect (with field- 
induced oxygen vacancies S{B), see Fig. 2) for two values 
of field-free deficiency parameter: S{0) ~ (solid line), 
and 5{0) =0.1 (dashed line). 

this implicit equation for the resulting magnetization 
m/ = M{B,d{B))/Ma{Q) is shown in Fig. | for the 
two values of zero-field deficiency parameter S{0). As is 
clearly seen, m/ exhibits a field- induced "fishtail" -like 
behavior typical for underdoped crystals with intragrain 
granularity (for symmetry and better visual effect we 
also plotted —mj in the same figure). The extra ex- 
tremum of the magnetization appears when the applied 
magnetic field B matches an intrinsic chemomagnetic 
field B^{5{B)) (which now also depends on B via the 
above-discussed magnetoconcentration effect). Notice 
that a " fishtail" structure of mf manifests itself even at 
zero values of field-free deficiency parameter (5(0) (solid 
line in Fig. ^) thus confirming a field-induced nature of 
intrinsic granularity ||, |l3, 14, 1^. At the same time, 
even a rather small deviation from the zero-field stoi- 
chiometry (with (5(0) =0.1) immediately brings about 
a paramagnetic Meissner effect at low magnetic fields. 
Thus, the present model predicts appearance of two in- 



terrelated phenomena, Meissner paramagnetism at low 
fields and "fishtail" anomaly at high fields. It would be 
very interesting to verify these predictions experimen- 
tally in non-stoichiomctric superconductors with pro- 
nounced networks of planar defects. 
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the subject. 
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